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Abstract 


Quantization  of  data  and  representa- 
^"rtion  of  the  coefficients  is  studied  in 
^^vdigital  matched  filters  for  weak-signal 
''^‘‘detection.  An  algorithm  for  optimum  co- 
£)^ef ficients  and  equations  for  the  optimum 
input  quantizer  are  obtained  for  the 
known  signal  in  additive  noise  problem. 
^^Some  numerical  performance  results  are 


given. 


I.  INTRODUCTION 


A matched  filter  is  often  used  as  a 
detector  for  testing  a hypothesis  about 
an  input  signal;  it  is  the  optimal  pro- 
cessor in  this  role  for  Gaussian  input 
noise,  and  may  also  be  considered  as  the 
optimum  processing  scheme  in  non-Gaussian 
noise  in  the  weak-signal  case.  We  will 
elaborate  on  this  in  the  next  section; 
the  general  theory  and  applications  of 
matched  filters  may  be  found,  for  exam- 
ple, in  [1,2]. 

In  this  paper  we  will  be  concerned 
with  digital  matched  filters  operating 
on  discrete- time  data,  and  we  will  exam- 
ine the  effects  of,  and  optimization  with 
respect  to,  finite-bit  representation  or 
quantization  of  the  coefficients  and  ana- 
log inputs.  Most  previous  investigations 
in  this  direction  have  been  concerned  ei- 
ther with  data  quantization  only  (e.g. 
[3])  or  have  assumed  only  simple  one-bit 
coefficient  representations  [2]. 

In  the  next  section  we  briefly  con- 
sider the  basic  results  on  local,  or 
weak-signal,  detection  of  signals  based 
on  quantized  data  and  finite-bit  coeffi- 
cients. In  Section  III  implementation 
of  optimum  digital  matched  filtering  is 
considered,  and  an  algorithm  for  deter- 
^ mining  the  best  coefficient  representa- 
2^“tion  is  discussed.  In  Section  IV  we 
give  some  numerical  examples. 

o 

C_}  II.  PERFORMANCE  WITH  QUANTIZED 

COEFFICIENTS  AND  INPUTS 


Let  us  assume  that  an  observation 
ZZZ vector  X = (X. ,X,, . . . ,X_)  is  available, 

i i — 1 z n 

and  is  described  by  the  equation 


0s,  + N, 


i-1,2,. 


,n,  0>O.  (1) 


Here  the  vector  s = (s^,s2, . . . ,sn>  is  a 

known  signal  vector,  9 is  the  amplitude 
of  the  signal,  and  the  vector 
N = (N1,N2, • . . ,NR)  is  a vector  of  inde- 
pendent identically  distributed  noise 
samples  each  with  symmetric  density  f. 

If  we  consider  as  a detection  sta- 
tistic for  testing  HQ:e=0  vs.  H^:9>0  a 


T = I g, (X.), 
i=l  1 1 


(2) 


and  use  as  a criterion  of  optimality  the 
differential  signal-to-noise  ratio 


DSNR  = 


[te^'e-o] 

Var{  T 1 , 


'6=0 


(3) 


we  find  (from  the  Schwarz  inequality) 
that  the  optimum  T maximizing  DSNR  is 


Topt  ~ gopt(Xi) 


where 


W<V  - -f  «i)/f(xi) 


(4) 


(5) 


the  prime  denoting  the  first  derivative 
of  the  function.  The  criterion  of  (3) 
is  a reasonable  one  for  the  weak-signal 
case,  being  a modified  case  of  the  usual 
SNR  criterion  [2].  It  is  also  well- 

known,  of  course,  that  T . is  the 

opt 

local ly-optimum  statistic  maximizing  the 
slope  of  the  power  function  for  testing 
Hq  vs. 


Hx  [4]. 

In  the  case  of  Gaussian  noise  T 


opt 

is  the  output  of  a linear  filter  matched 
to  the  signal  vector  s.  In  the  general 
case,  we  may  consider  T0pt  to  be  the  out- 
put of  a similar  filter  preceded  by  the 


instantaneous  nonlinearity  g 


opt  * 


For  the 


digital  matched  filter  both  the  coeffi- 
cients, the  s^,  and  the  data  inputs,  the 

gQpt(X^),  have  to  be  replaced  by  suitably 

quantized  versions  r.  and  q(X-),  respec- 
tively. 


' ease; 
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To  optimize  the  performance  of  the 
digital  matched  filter,  then,  we  have 
to  optimize  DSNR  of  (3)  for  the  case 
where  the  detection  statistic  is 


T = E r.  q(X.),  (6) 

d i=l  1 1 

the  optimization  being  with  respect  to  a 
choice  of  the  r,  and  the  input  quantizer 

u * 

q,  given  2 =2k  levels  or  b bits  for  co- 
efficient representation  and  2^=2m  lev- 
els or  t bits  for  the  data  quantization. 
Note  that  (6)  represents  a general  dig- 
ital matched  filter  operating  on  quan- 
tized data,  for  any  kind  of  input  noise 
density. 

Applying  (3)  to  T^,  we  find  that 
DSNR  for  with  an  odd-symmetric, 
even- state  quantizer  q is  given  by 


DSNR.  = 
a 


( 


2 [ij 


E i.z[P(a.  . ) -F (a  • ) ] 


where  F is  the  distribution  function 
corresponding  to  f.  In  (7)  the  are 

the  output  levels  for  positive  x of  the 
symmetric  quantizer  q,  with  q(x)  being 
tj  whenever  x is  in  the  interval 

(aj,  »j_^) • The  breakpoints  a^  satisfy 

a.  < a.  ,,  and  by  definition  a_=0  and 
3 j-1  n» 

a0-». 

The  problem  is  to  maximize  the 
quantity 

- n ,2 


U Vi] 


which  is  a factor  in  DSNR. 


In  addition. 


for  a given  set  of  levels  t...  for  the 

quantizer  q,  the  best  set  of  breakpoints 
Aj  may  be  obtained.  Then  the  overall 

performance  of  the  digital  matched  filter 
may  be  optimized  with  respect  to  alloca- 


tion of  b+t  bits  between  coefficient  and 
data  quantization. 

III.  PERFORMANCE  OPTIMIZATION 

As  discussed  in  the  previous  section, 
the  objective  is  to  maximize  DSNR^  of  (7). 

The  two  factors  in  (7)  are  decoupled,  so 
that  the  coefficient  representation  and 
data  quantization  problems  can  be  treated 
independently  of  each  other.  We  first 
consider  maximization  of  the  coefficient 
factor  J defined  in  (8) . This  quantity 
can  be  written  as 

(r-s)2 

J - (9) 


where  s is  the  previously  defined  signal 

vector  and  r = (r, ,r.,...,r  ),  the  vector 
— l z n 

of  filter  coefficients  by  which  the  ref- 
erence signal  s is  represented.  If  we 
let  <f>  be  the  angle  between  s and  r,  so 
that 


l|r||2||s||2 
we  can  express  J as 

J = | | s | | 2cos2<}>  . 


Thus,  maximization  of  J for  a given  ref- 
erence signal  is  simply  maximization  of 
2 

cos  <t>,  and  we  therefore  have  to  pick  the 
vector  r closest  to  s. 

Since  it  was  not  easy  to  aet  an 
analytic  solution  for  the  coefficient 

2 

vector  r maximizing  cos  <f>  for  a given  s 
(and  with  a constraint  on  the  number  of 
bits,  b,  for  coefficient  representation), 
an  efficient  computer  technique  was  de- 
veloped which  is  described  below. 

It  was  assumed  that  the  coefficient 
representation  r is  obtained  through  a 
quantizer  with  one  of  three  possible 

ranges  of  2b=2k  levels: 

(a)  The  levels  0,l,...,2k-l 

(b)  The  levels  -k+1,  -k+2,...,k 

(c)  The  levels  -k,  -k+l,...,k-l 

Range  (a)  is  obviously  to  be  used  if  s 
has  positive  components.  A signal  vector 
with  all  negative  components  may  be  com- 
plemented, and  therefore  be  represented 
in  this  range  also.  Ranges  (b)  and  (c) 
are  more  natural  choices  if  positive  as 


well  as  negative  components  are  present 
in  s.  Note  that  we  are  considering  the 
class  of  quantizers  with  an  even  number 
of  levels,  one  of  which  is  the  0 level. 

The  search  algorithm  for  r is  des- 
cribed by  the  following  sequence  of 
steps.  It  is  assumed  that  the  signal 
vector  components  have  been  ordered, 
that  is  s^<s2£. . .<sn,  and  are  non-zero. 

0.  Initialization:  j * 2k-l,  00. 

1.  Check  if  sn>0.  If  not,  invert  the 

sequence  of  components  in  £ so  that 
now  s = (-s  ,-sn_1,...,-SjT»  and  9° 

to  step  3. 

n 

2.  Check  if  I s->0.  If  the  sum  is 

i=l  1_ 

negative,  invert  the  sequence  of 

components  in  s,  so  that  now 

s = (-s  ,-s  .,...,-3.).  Note  that 

— n n- 1 X 

inverting  the  sequence  in  step  1 

results  in  a sequence  which  fulfills 

the  condition  in  this  step. 

3.  Form  the  vector  s(j)  = ( s -^  ( j ) , 
s2(j) , . . . ,sn(j) ) where  si ( j ) =jsi/sn. 
Let  I^(j)  be  the  largest  integer 
less  than,  or  equal  to,  s^( j).  Thus 
In(j)  = j;  the  algorithm  now  looks 

for  the  first  (n-1)  components  of 
r ( j ) = (rx ( j ) ,r2 (j) , . . . , rn_1 ( j ) , j) 

so  that  the  cosine  squared  of  the 
angle  between  s and  r(j)  is  maxi- 
mized. 

4.  If  j=k,  go  to  step  6.  If  j<k,  go 
to  step  7. 

5.  If  s^(j)>0,  consider  the  two  possi- 
bilities ( j ) and  U^(j)  for  r^tj), 
where  IT  (j)  is  the  smallest  integer 
larger  than,  or  equal  to,  si(j). 

If  si(j)<0,  set  ri(j)=0.  Go  to 
step  8 . 

6.  If  s^(j)  > -k+1,  consider  the  pos- 
sibilities I^(j)  and  ( j ) for 
ri(j).  If  s±  ( j ) < k+1,  set 

r^(j)  * -k+1.  Go  to  step  8. 

7.  If  s^(j)  £ -k,  consider  the  possi- 
bilities I^fj)  and  ( j ) for  r ^ ( j ) . 
If  ( j ) < -k,  set  r±(j)  * -k. 


8.  Compute  the  square  of  the  cosine, 

2 

c (j),  of  the  angle  between  s and 
r(j)  for  each  different  possible 
combination  of  components  of  r(j); 
let  r (j)  be  the  vector  yielding 

—Til  2 2 

a maximum  value  c (j)  for  c (j). 

2 m 2 
If  c_  > C,  assign  the  value  c_(j) 
m in 

to  C and  define  r to  be  r (j)  . 

— — m 

9.  If  j>l,  set  j=j-l  and  go  to  step  3. 

10.  Stop.  The  optimum  coefficient  vec- 
tor is  r,  and  C is  the  square  of 
the  cosine  of  the  angle  between  r 
and  £,  if  the  two  conditions  of 
steps  1 and  2 had  been  fulfilled. 
Otherwise  the  optimum  coefficient 
vector  is  in  inverse  order  in  r. 

On  an  intuitive  basis,  one  might 
choose  the  best  coefficient  vector  in 
the  following  way  for  the  simple  case 
where  the  s^  are  positive.  Pick  the 

vector  of  coefficients  r,  where  the  com- 
ponents r^  are  the  integers  closest  to 

the  s^(2k-l).  However,  it  is  clear  that 

r need  not  be  an  optimum  set  maximizing 
2 

cos  (J>.  This  can  be  seen  by  considering 
the  case  b=l,  giving  two-level  represen- 
tation, for  (n-1)  identical,  small  s^ 

and  a large  value  for  sn>  in  the  next 

section  this  is  further  illustrated  for 
an  example  which  is  not  as  extreme. 


We  also  need  to  consider  optimiza- 
tion with  respect  to  the  data  quantizer. 
For  the  symmetric  even-state  input  quan- 
tizer we  are  considering,  specification 

of  the  number  of  levels  2t=2m  (with  t 

bits)  fixes  the  levels  at  values 

+l,+2, . . . ,+m.  With  =m-j+l,  j=l,2, . . . ,m, 

the  second  factor  K in  (9) , defined  by 


m 


Z l/[F(a.  . ) -F  (a  . ) ] 
1=1  3 J"1  3 


(12) 


can  be  maximized  with  respect  to  the  a ^ . 

The  optimum  set  can  easily  be 

obtained  for  specific  densities  f.  Set- 
ting the  partial  derivative  of  K with 
respect  to  a^  equal  to  zero,  we  find  a 


necessary  condition  for  the  maximizing 
values  of  a j : 


4 


i J 
1 


m 


-f  (a,) 
TTIj 


■E 


+ yn  ^Aj[f  (aj>-f  (aj-i)] 


1 m 

I A.“[F(a.  . ) -F (a .)] 
j=l  3 31  3 

3*1# • • • i®*l 
(13) 


Equation  (13)  may  be  solved  for  specific 
densities.  For  the  case  of  a Gaussian 

2 

density  with  variance  a , the  equations 
reduce  to 


3 r>ii 


aj 


L,  j = l. 


where 


L 


l A-Jf  (a.)-f  (a.  ,)  ] 
1=1  3 3 31 
m . 

E A.z[F(a.  . ) -F (a . ) ] 
j=l  3 3-1  3 


(14) 


(15) 


Substitution  of  (14)  in  (15)  leads  to  a 
single  equation  for  L,  which  may  then  be 
solved  (numerically)  and  hence  the  opti- 
mum a^  can  be  obtained.  The  next  sec- 
tion contains  some  specific  numerical 
results. 


IV.  RESULTS  AND  DISCUSSION 

The  considerations  of  the  previous 
section  were  applied  to  several  specific 
cases,  two  of  which  are  presented  here. 

A signal  vector  s of  10  components  and 
one  of  15  components  is  shown  in  Table  I, 
together  with  optimum  representations 
for  one,  two  and  three  bits.  The  maxi- 
mum values  J of  J [Equation  (11) ] 

max  ~ 

normalized  by  | |s|  | are  also  shown  in 
Table  I.  The  results  indicate  that  even 
with  only  two  bits  for  coefficient  rep- 
2 

resentation  (cos  #)  is  very  close  to 
max 

unity,  which  corresponds  to  the  analog 
case.  The  exact  numerical  values  depend 
on  the  coefficient  vector  and  its  length. 

In  the  second  part  of  Section  III 
the  optimization  of  the  second  factor  K 
in  (7)  was  discussed.  The  optimum  per- 
formance for  the  specific  case  of 
Gaussian  noise  with  unit  variance  is 
shown  in  Table  II.  Again,  in  this  case 
the  value  K=1  is  obtained  for  unquan- 
tized data.  Table  III  combines  the  re- 
sults of  Tables  I and  II,  giving  the 
overall  performance  for  the  different 
allocations  of  a fixed  number  of  bits 
(six  bits)  between  coefficient  represen- 


tation and  data  quantization.  In  general, 
the  optimum  allocation  will  depend  on  the 
length  of  the  vector  s,  its  components 
and  the  type  of  noise  density.  It  is 
seen  that  one  need  consider  only  a small 
number  of  bits  to  achieve  near-analog 
performance. 

In  Section  III  an  approximation 
method  of  quantizing  the  coefficient 
sequence  s was  also  discussed.  Table  IV 
gives  an  example  of  this  approximation 
method  using  the  same  sequence  as  in 
Table  1(a)  with  n=10.  Note  that  the  dis- 
crepancy is  larger  for  the  lower-order 
quantizers,  as  expected.  In  fact,  the 
result  of  the  approximation  for  3 bits 
is  identical  to  that  in  Table  1.  It  is 
seen  that  in  this  example,  the  approxi- 
mation methods  gives  performance  very 
close  to  optimum. 
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Table  I:  Optimum  Coefficient  Representation 

(a)  n=10 


COEFFICIENTS 

-1.82  -0.92  -0.08  0.77  1.22 

1.58  1.98  2.30  2.41  3.10 

J 

max  , 2.x 

5-  = (cos  4>) 

||s||2 

ONE  BIT 

REPRESENTATION 

0 0 0 0 1 
11  1 1 1 

0.792 

TWO  BIT 

REPRESENTATION 

-1  -1  0 11 

1 2 2 2 2 

0.957 

THREE  BIT 
REPRESENTATION 

-2-1  0 12 
2 3 3 3 4 

0.990 

(b)  n=15 

COEFFICIENTS 

J 

max  . 2 . » 

=-  = (cos  $) 

l|s||2 

ONE  BIT 

REPRESENTATION 

0 0 0 0 0 
0 0 0 0 0 
11111 

0.604 

TWO  BIT 

REPRESENTATION 

-1  -1  -1  -1  -1 

0 0 0 1 1 

1 2 2 2 2 

0.907 

THREE  BIT 
REPRESENTATION 

-3  -1  -1  -1  -1 

0 0 0 1 1 

1 2 2 3 3 

0.979 
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